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' Abstract 

^ ' We review the hypergeometric function approach to Feynman diagrams. Special con- 

, sideration is given to the construction of the Laurent expansion. As an illustration, we 
describe a collection of physically important one-loop vertex diagrams for which this ap- 

, proach is useful. 

r~| ' 1. Introduction. Recent interest in the mathematical structure of Feynman diagrams 

"V* ' has been inspired by the persistently increasing accuracy of high-energy experiments and the 

advent of the LHC epoch. For stable numerical evaluation of diagrams, a knowledge of their 
r~| I analytical properties is necessary. We will review some of the progress in this area, focusing on 

the hypergeometric function representation of Feynman diagrams. 
Ti^^j- . Forty-five years ago, Regge proposed [T] that any Feynman diagram can be understood 

^ I in terms of a special class of hypergeometric functions satisfying some system of differential 

^ equations so that the singularity surface of the relevant hypergeometric function coincides 

^SJ ■ with the surface of the Landau singularities [2] of the original Feynman diagram0 Based on 

I Regge's conjecture, explicit systems of differential equations for particular types of diagrams 

' have been constructed. For some examples, the hypergeometric representation for A'^-point 

. one-loop diagrams has been derived in Ref. [1] via a series representation (Appell functions 

i and Lauricella functions appear here), the system of differential equations and its solution in 

^ ' terms of Lappo-Danilevsky functions \5] has been constructed in Ref. [6], and the monodromy 

structure of some Feynman diagrams has been studied in Ref. [7] . 

A review of results derived up to the mid-1970's can be found in Ref. It was known at 
I that time that each Feynman diagram is a function of the "Nilsson class." This means that the 

Feynman diagram is a multivalued analytical function in complex projective space CP". The 
singularities of this function are described by Landau's equation. Later, Kashiwara and Kawai 
showed [9] that any regularized Feynman integral satisfies some holonomic system of linear 
differential equations whose characteristic variety is confined to the extended Landau variety. 

The modern technology for evaluating Feynman diagrams is based mainly on techniques 
which do not explicitly use properties of hypergeometric functions, but are based on relation- 
ships among the Feynman diagrams derived from their internal structure!! It was shown, for 

*e-mail: kalmykov@theor.jinr.ru, kniehl@mail.desy.de, BFL_Ward@baylor.edu, scott.yost@citadel.edu 
^ For a review of different approaches to the analysis of the singularities of Feynman diagrams see Ref. ^ . 
^By "internal structure," we mean any representation described in standard textbooks, such as Ref. [10) . 
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example, that there are algebraic relations between dimensionally regularized [II] Feynman 
diagrams with different powers of propagator [12j. Tarasov showed in 1996 that similar al- 
gebraic relations could also be found relating different dimensions of the integral [H]. The 
Davydychev- Tarasov algorithm [131 [H] allows a Feynman diagram with arbitrary numerator 
to be transformed into a linear combination of diagrams of the original type with shifted pow- 
ers of propagators and space-time dimension, multiplied by a linear combination of tensors 
constructed from the metric tensor and external momenta. This set of algebraic relations is 
analogous to contiguous relations for hypergeometric functions]! 

Solving the algebraic relations among Feynman diagrams allows them to be expressed in 
terms of a restricted set called "master integrals." Such a solution is completely equivalent to 
the differential reduction of hypergeometric functions [16[ [T71 [T8] . The technique of describing 
Feynman diagrams by a system of differential equations was further extended in Ref. [19j , where 
it was realized that the solution of the recurrence relations can be used to close the system of 
differential equations for any Feynman diagram. This led to useful techniques for evaluating 
diagrams [201 [2T] . Most of the progress to date in this type of analysis has been for diagrams 
related to the "Fuchs" type of differential equation, with three regular singular points [22j^ . 




Since Feynman diagrams are often UV- or IR-divergent, it is important to also consider 
the construction of the Laurent expansion of dimensionally-regularized diagrams about integral 
values of the dimension (typically d = 4 — 2e:). This is called an "e expansion" of the diagram. 
For practical applications, we need the numerical values of the coefficients of this expansion. 
Purely numerical approaches are under development {e.g. Ref. [21]), but this is a complex 
problem for many realistic diagrams having UV and IR singularities and several mass scales. 

The case of one-loop Feynman diagrams has been studied the most. The hypergeometric 
representations for N-point one-loop diagrams with arbitrary powers of propagators and an 
arbitrary space-time dimension have been derived for non-exceptional kinematic^ by Davy- 
dychev in 1991 [25j. His approach is based on the Mellin-Barnes technique |26[ . The results 
are expressible in terms of hypergeometric functions with one less variable than the number of 
kinematic invariants. 

An alternative hypergeometric representation for one-loop diagrams has been derived re- 
cently in Ref. [28], using a difference equation in the space-time dimension. This approach has 
been applied only to a set of master integral^, but, fortunately, an arbitrary A^-point function 
can be reduced to the set of master integrals analytically [29l[30]. In Ref. [28j, the one-loop 
A^-point function was shown to be expressible in terms of hypergeometric functions of A^ — 1 
variables. One remarkable feature of the derived results is a one-to-one correspondence between 
arguments of the hypergeometric functions and Gram and Cayley determinants, which are two 
of the main characteristics of diagrams. 

Beyond one loop, a general hypergeometric representation is available only for sunset-type 
diagrams with arbitrary kinematics [31], with a simpler representation for particular kinematics 
[S21E3]. In all other cases beyond one loop, master integrals have been expressed in terms of 
hypergeometric functions of type p-Fp-i fMj- 

The program of constructing the analytical coefficients of the e-expansion is a more compli- 
cated matter. The finite parts of one-loop diagrams in d = 4 dimension are expressible in terms 

''The full set of contiguous relations for generalized hypergeometric functions pFq is found in Ref. [15| . 
*The analysis of some diagrams with four regular singularities was done recently in Ref. |23j. 
^"Non-exceptional kinematics" refers to the case where all masses and momenta are non-zero and not pro- 
portional to each other. 

^The hypergeometric representations of one-loop master integrals of propagator and vertex type have been 
constructed in ,26. i27] . 
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of the Spence dilogarithm function [35J. However, only partial results for higher-order terms 
in the e-expansion are known at one loop. The all-order e-expansion of the one-loop propa- 
gator with an arbitrary values of masses and external momentum has been constructed [3_7J 
in terms of Nielsen polylogarithms [36j . The term linear in e for the one-loop vertex diagram 
with non-exceptional kinematics has also been constructed in terms of Nielsen polylogarithms 
|38j . It was shown in Ref. [39] that the all-order e expansion for the one- loop vertex with 
non-exceptional kinematics is expressible in terms of multiple polylogarithms of two variables 

m- 

Beyond these examples, the situation is less complete. The term linear in e for the box 
diagram is still under construction. Some cases for particular massed have been analyzed 
[41^ I42j. Many physically interesting particular cases have been considered beyond one loop. 
Among these are the e expansion of massless propagator diagrams and the sunset diagram 



2. Hypergeometric Functions. Let us recall that there are several different ways to de- 
scribe special functions: (i) as an integral of the Euler or Mellin-Barnes type; (ii) by a series 
whose coefficients satisfy certain recurrence relations; (iii) as a solution of a system of differ- 
ential and/or difference equations (holonomic approach). These approaches and interrelations 
between them have been discussed in series of a papers [35]. In this section, we review some 
essential definitions relevant for each of these representations. 

• Integral representation: An Euler integral has the form 

^{aJ,P)= I UiP,{xu--- ,Xkf^x'^' ■■■x'^''dxi---dxk , (1) 
Js 

where Pi is some Laurent polynomial with respect to variables xi, • • • , x^: Pi{xi, • • • , x^) = 
J2 Cuji-uikXi^ ■ ■ ■ ^ with S Z, and Oj, /3j £ C. We assume that the region S is chosen 
such that the integral exists. 

A Mellin-Barnes integral has the form 

<!> (Ojs, Cj, dj, 7, x) = / dz\ . . . dZra^^ -p i^m T \ TT"^! ^ • • • -^m ) (^) 

J7+iR -'--'-fc=l-'- \1^T=\ '^krZr + "fcj 

where ajs,bkr,Ci,dj € M, € C, and 7 is chosen such that the integral exists. This 
integral can be expressed in terms of a sum of the residues of the integrated expression. 

• Series representation: We will take the Horn definition of the series representation. 
In accordance with this definition, a formal (Laurent) power series in r variables, 

$(f) = ^C'(m)x^ = C(mi,m2,--- ,m^)x™i (3) 

m,i,m2,--- ,mr 

is called hypergeometric if for each i = 1, • • • , r the ratio C{rh+ei)/C{rh) is a rational func- 
tiorll in the index of summation: rh = (mi, • • • , rrir), where e*,- = (0, • • • , 0, 1, 0, • • • , 0), is 
unit vector with unity in the j^^ place. Ore and Sato [46J found that the coefficients of 
such a series have the general form 

C(m) = nLiAri?(m) f nf=ir(M,(m) + 7i + 1)^ , (4) 



^In Ref. [28j . box diagrams have been written in terms of the Lauricella-Saran function Fs of three variables, 
and a one-fold integral representation was established for their all-order e expansion. However, it is not proven 
that this representation coincides with multiple polylogarithms. 

*A "rational function" is any function which can be written as the ratio of two polynomial functions. 
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where iV > 0, Aj,7j G C are arbitrary complex numbers, fij : V ^ 'L are arbitrary 
linear maps, and R is an arbitrary rational function. The fact that all the F factors are 
in the denominator is inessential: using the relation r(z)r(l — z) = vr/ sin(7rz), they can 
be converted to factors in the numerator. A series of this type is called a "Horn-type" 
hyper geometric series. In this case, the system of differential equations has the form 

^fcT^ I Z-^(^) ( >J^fc£- 1 ^(^) ' j = l,---,r, (5) 



(6) 



where Pj and Qr are polynomials satisfying 




C{m + ej) _ Pj{m) 



C(m) Qj{7ri) 



ti satisfies the ec 
j = l,---,r, (8) 




• Holonomic representation: A combination of differential and difference equations can 
be found to describe functions of the form 

'^(^>^'^)= t fn^^-T-rW-wIirln,^^,^, (7) 

ki,--- ,kr=0 

where 14^ is an r x m matrix. In particular, this function satisfies the equations 
d^{z,x,W) 

^ |^Zi$ + ^iy,x,^ I =0, i = l,---,m, (9) 
where ujj is the j*^ column of the matrix W . 

3. Construction of the all-order e expansion of hypergeometric functions. Recently, 
several theorems have been proven on the all-order e expansion of hypergeometric functions 
about integer and/or rational values of parameters [321 EH HZl HHl 113 EOl ED E^- For hyper- 
geometric functions of one variable, all three of the representations (i)-(iii) described in the 
previous section are equivalent, but some properties of the function may be more evident in 
one representation than another. 

In the Euler integral representation, the most important results are related to the construc- 
tion of the all-order e expansion of Gauss hypergeometric function with special values of pa- 
rameters in terms of Nielsen poly logarithms [37j. There are several important master integrals 
expressible in terms of this type of hypergeometric function, including one-loop propagator- 
type diagrams with arbitrary values of mass and momentum [26], two-loop bubble diagrams 
with arbitrary values of masses, and one-loop massless vertex diagrams with three non-zero 
external momenta |53j . 

The series representation (ii) is an intensively studied approach. The first results of this 
type were derived in the context of the so-called "single-scale" diagrams [54J related to multiple 
harmonic sums. These results have been extended to the case of multiple (inverse) binomial 
sums [STj that correspond to the e-expansion of hypergeometric functions with one unbalanced 
half-integer parameter and values of argument equal to 1/4, or diagrams with two massive- 
particle cuts. Particularly impressive results involving series representations were derived in 
the framework of the nested-sum approach for hypergeometric functions with a balanced set 
of parameters in Refs. [471 148j . [f| and in framework of the generating- function approach for 
hypergeometric functions with one unbalanced set of parameters in Refs. [33 1 l5T | [58 | [59]. 



Computer realizations of nested sums approach to expansion of hypergeometric functions are given in | 



4 



An approach using the iterated solution of differential equations has been explored in Refs. 
[33| SHI [Sni [52] . One of the advantages of the iterated-solution approach over the series approach 
is that it provides a more efficient way to calculate each order of the e expansion, since it relates 
each new term to the previously derived terms, rather than having to work with an increasingly 
large collection of independent sums at each order. This technique includes two steps: (i) 
the differential-reduction algorithm (to reduce a generalized hypergeometric function to basic 
functions); (ii) iterative solution of the proper differential equation for the basic functions 
(equivalent to iterative algorithms for calculating the analytical coefficients of the e expansion) . 

An important tool for constructing the iterative solution is the iterated integral defined as 
I{z; ak,ak-i, . . . , ai) = Jq j^^I{t; a/c-i, . . . , ai) , where we assume that all aj 7^ 0. A special 
case of this integral, 

Gmk,m,,_l,...,m,x{z]ak-, ■ ■ ■ ,ai) = I{z; 0, • - ; )flfc) 0; • ^- ' Q jflfc-l)""") 0) "J ''^^ '"^l); 

times mfc_i-i times mi-i times 

where all 7^ 0, is related to the multiple polylogarithm [301 EH 

Lifci,fc2,...,fc„(xi,X2,...,a;„) = ^^^•••^^ (10) 

m„>m„_i>- >m2>mi>0 "^1 "^2 

by 



'-^m„,m„_i,...,mi [Z'j Xm Xn—lj ■ ■ ■ , Xi) — I) i-'lmi,m2,.--,"in-l,"^n I ; i ■ ■ ■ i j 

\Xi X2 Xfi—l Xn ^ 

Lifei,fc2,.-,fcn-i,fcn (yi: y2, • • • , yn-i,yn) = (-i)"'Gfc„,fc„_i,...,fc2,fei 1; — , • • • , — — 

V Vn yn X • • • X 2/1 

In Eq. ([To]), k = ki + k2 + ■ ■ ■ + kn is cahed the "weight" and n the "depth." Multiple 
polylogarithms pO|) are defined for < 1 and \zi\ < l{i = 1, . • • • ,n — 1) and for |z„| < 1 if 
m„ < 2. We mention also that multiple polylogarithms form two Hopf algebras. One is related 
to the integral representation, and the other one to the series. 

A particular case of the multiple polylogarithm is the "generalized polylogarithm" defined 

by 



• • - JL 

— — -r ^ =Lifcj,fc2,...,fc„ (1,1,--- ,1,2:) , (11) 

mn>mn-\>--->m\>0 L 2 " 

where j^l < 1 when all A;^ > 1, or |z| < 1 when fe„ < 2. Another particular case is a "multiple 
polylogarithm of a square root of unity," defined as 

V si,S2, ■>sn / ^ ^ ^ r, '"''^ '"'1 

^ m„>m„-i>--mi>0 ^ 

where s = (si,---Sn) and a = {ai,--- ,an) are multi-indices and at belongs to the set of 
the square roots of unity, = ±1- This particular case of multiple polylogarithms has been 
analyzed in detail by Remiddi and Vermaseren [62 J^ . 

Special consideration is necessary when the last few arguments afc_j,afc_j_i, ™- the 

integral I{z; ai, • • • , a^) are equal to zero, which is called the "trailing-zero" case. It is possible 
to factorize such a function into a product of a power of a logarithm and a multiple polylog- 
arithm. An appropriate procedure for multiple polylogarithms of a square root of unity was 



^"As was pointed out by Goncharov [40], the iterated integral as a function of the variable z has been studied 
by Kummer, Poincare, and Lappo-Danilevky, and was called a hyperlogarithm. Goncharov [ID] analyzed it as a 
multivalued analytical function of ai, . . . , a^, 2. From this point of view, only the functions considered in Ref. 
[63] are multiple polylogarithms of two variables. 
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described in Ref. [b2] and extended to the case of multiple polylogarithms in Ref. [63]. For the 
numerical evaluation of multiple polylogarithms or its particular cases, see Ref. [MIES]. Let 
us consider the Laurent expansion of a generalized hypergeometric functions of one variable 
pFp-i{A; B; z) with respect to its parameters. Such an expansion can be written as 

pFp^i{A; B; z) = pFp_i(Ao; Bq; z) 



mi,l 



= pFp.i{Ao;Bo;z)+ ^ UP^,UPrl(Ai - Aoi^iB, - B^jY^ L^^g{z) , (13) 

mi,lj=l 

where pFp-i{A; B; z) is a hypergeometric function defined by pFp^i{A; B; z) = YlT=o np-i/p'l^ fr 

and {A)j is the Pochhammer symbol: {A)j = T(A + j)/T{A). Our goal is to completely de- 
scribe the coefficients g{z) entering the r.h.s. of Eq. p^ . To reach this goal, we must first 
characterize the complete set of parameters for which known special functions suffice to express 
the coefficients. Beyond this, we wish to identity the complete set of new functions which must 
be invented in order to express all of the coefficients in the Laurent expansion. 

The first simplification comes from the well-known fact that any hypergeometric function 
pFp-i{A + m; B + k; z) may be expressed in terms of p other functions of the same type: 

p 

Rp+i{A,B,z)pFp.i{A + m;B + k;z) = ^Rj{A,B,z)pFp^i{A + ek;B + Ek;z) , (14) 

i=i 

where rh, k, Cfc, and Ek are lists of integers, and the Rk are polynomials in the parameters A, B, 
and z. In particular, we can take the function and its first p—1 derivatives as a basis for the 
reduction (see Ref. [16] for the details of this approach). Then Eq. ()14p will take the fornix 

~ ^ ~ / d \ 

Rp+i{A,B,z)pFp_iiA + m;B + k]z) = J2Rk{A,B,z)y—j pFp_i{A; B; z) , (15) 

with a new polynomial In this way, the problem of finding the Laurent expansion of the 
original hypergeometric function is reduced to the analysis of a set of basic functions and the 
Laurent expansion of a (formally) known polynomial. 

As is well known, hypergeometric functions satisfy the differential equatiorEl 



pF,_i(A;B;l) = 0. (16) 



Due to the analyticity of the hypergeometric function p-Fp_i {A; B; z) with respect to its parame- 
ters Ai,Bk, the differential equation for the coefficients ^{z) of the Laurent expansion could 
be directly derived from Eq. p6|) without any reference to the series or integral representation. 
This was the main idea of the approach developed in Refs. [331 H9| [50l [52| I60j. An analysis of 
this system of equations and / or their explicit analytical solution gives us the analytical form of 
g{z). It is convenient to introduce a new parametrization, Ai A^i + Oje, Bj — > B^i + biS , 
where e is some small number, so that the Laurent expansion ()13p takes the form of an "e 
expansion," 

oo oo 

pFp_i(l+ae;^ + 6e;z) =pFp_i(i';B;^) + ^e'=L.g_^(z) = ^e%^g^^(z) , 

k=l k=0 



^^For simplicity, we will assume that no diflFerence Bk — Aj is a positive integer. 

This equation follows from Eqs. @ - ®, where P{j) = n^^-^(Afc + j) and Q{j) = (j + l)n^^}(Bfc + j) 
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where L^j^^^{z) = pFp-i(A; B; z). The differential operator can also be expanded in powers of 
e: 



where 6 
and 
i(p) . 



1 



nf^i {0+Ai+aie)--eUlz\ {e+Bk-l+bke) 



p 



j=0 



iA,B,a,b,z), (17) 



Zj^ , the upper index gives the order of the differential operator, D^^ = n^^-^Ofc 







^+Ai)--eui-\{e+Bk-i) 
p-i 



(18) 



k=l 



X,{A,B)--Y,{A,B) 



QP-^-j 



where Xj{A, B) and YjiA, B) are polynomials. Combining all of the expansions together, we ob- 



tain a system of equations J^^g X^^o ^^j' ''^ ^dbr- •(-^) ~ ^ ' which could be split into follow- 



ing system (each order of e): (e") D^^^L^^^iz) = ; {e", 1 < k < p) J2r=o "'L^ ^ 



d,b,r—j * 



(p-k) 



k—r ^ 



^(p-k) 



. Further simplification comes from the explicit 



forms of and the polynomials Xj{A, B),Yj{A, B) in Eq. (I19p . For example, for integer 

values of parameters, we can put = 0, Bj = 1, so that all of the Xj{A, B) and YjiA, B) are 
equal to zero. Further details can be found in our pape rs, Refs. pSj [50l EH [60] . 
Here, we will mention some of the existing results. |^ 

• If Ii,l2,l3 are arbitrary integers, the Laurent expansions of the Gauss hypergeometric 
functions 



2Fi(/i + ae, h + 6e; Is + I + ce; z) 
2Fi(/i + £ + ae, /2 + he; h + ce; z) 



2Fi(/i + 2 + ae, I2 + 2 + 6e; Is + f + ce; z) 
2Fi(/i + I + ae, I2 + be; I3 + f + ce; z) 



are expressible in terms of multiple polylogarithms of arguments being powers of g-roots 
of unity and a new variable, that is an algebraic function of with coefficients that are 
ratios of polynomials. 

• \{ A^B are lists of integers and /, p, q are integers, the Laurent expansions of the gener- 
alized hypergeometric functions 

pFp_i (1 + ae, 2 + /; + te; z) , pFp.i{A + ae; B + be,^ + I; z) 

are expressible in terms of multiple polylogarithms of arguments that are powers of q- 
roots of unity and a new variable that is an algebraic function of z, with coefficients that 
are ratios of polynomials. 

• If A,B are lists of integers, the Laurent expansion of the generalized hypergeometric 
function 

pFp^i{A + ae; B + be; z) 
are expressible via generalized polylogarithms (jlip . 

We should also mention the following case [48J in which the e expansion has been constructed 
via the nested sum approach: 

If p, g, Ik are any integers and A, B are lists of integers, the generalized hypergeometric function 
pFp_i({2+i'+ae}^,/;+ce;{2+5+6e}^,/2+(fe;^) 



In the following, we will assume that a, b, c are an arbitrary numbers and e is a small parameter. 
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is expressible in terms of multiple polylogarithms of arguments that are powers of g-roots of 
unity and the new variable z^/'^, with coefficients that are ratios of polynomials. A hypergeo- 
metric function of this form is said to have a zero-balance set of parameters. 

We will now demonstrate some algebraic relations between functions generated by the e 
expansion of hyper geometric functions with special sets of parameters. Let us consider the 
analytic continuation of the generalized hypergeometric function p+iFp with respect to the 
transformation 2; — > 1/z j34j : 

P+i 

ai, 02, • • • , ap+i 



1 



bi,b2, • • • ,bp 



E 

k=l 

-b2, 



-cik) 



n^=ir(5,-afc) 



• • • , l + ak-bp 

l + Ofc — Op+i 



(19) 



1, Ofc, l+ak-bi, 1+afc- 
l + Ofc-ai, l + afc-a2, 
where none of the differences between pairs of parameters aj — ak is an integer. 

On the r.h.s. of Eq. ()19p . we actually have a hypergeometric function p+iFp, since one of 
the parameters is always equal to unity. If we make the replacements 



q 



r 

q 



- + bje 



in Eq. ()19p . we obtain the relation 



p 

E 



\ + &, {1 + aje}p 
{1 + 5,4 



(20) 



where the Cr are constants. Another relation follows if we choose in Eq. ()19p the following set 
of parameters: 



j = 1, - • • ,p + 1 , bk ^bkS , k = 1, 



1 



- +bpe . 



Then we have 



{aje} 



p~l ' q 



p+1 

+ b„s 



p 

E 

s=l 



1 - ^ - ce, {1 + aje}p 



(21) 



where the c are constants. In this way, we find a proof of the following result: 
Lemma: When none of the difference between two upper parameters is an integer, and the 
differences between any lower and upper parameters are positive integers, the coefficients of the 
£ expansion of the hypergeometric functions 

I _ / A+ae 



p+iFp 



A+^+ae 
B + '- + be 



B+b£,I+^+ce 



-iFp 



I+^+ce,A+ae 
B+be 



where A, B,a,b,c and I are all integers, are related to each other. 

Note that none of the functions of this lemma belongs to the zero-balance case. 

4. One-loop vertex as hypergeometric function. Let us consider now the one-loop 
vertex diagram. We recall that any one-loop vertex diagram with the arbitrary masses, external 
momenta and power of propagators can be reduced by recurrence relations to a vertex-type 
master integral (with all powers of propagators being equal to unity) or, in the case of zero 
Gram and/or Cayley determinants, to a linear combination of propagator- type diagrams [29j . 
In the case of non-zero Gram and/or Cayley determinants, the one- loop master integrals are 
expressible in terms of linear combinations of two Gauss hypergeometric functions and the 
Appell function Fi [271[28]. 

Our starting point is the hypergeometric representation for one-loop diagrams with three 
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c, 




C4 

Q 

^\ m 



m 







^\ 

J. 



C5 

Q 

m /\ m 









Qi 

/\ m 



m 



Figure 1: One- loop vertex- type diagrams expressible in terms of generalized hypergeometric functions. 
Bold and thin lines correspond to massive and massless propagators, respectively. 



arbitrary external momenta and one massive line or two or three massive lines with an equal 
masses, derived in Ref. [26]. Let us consider a one-loop vertex-type diagram, as shown in Fig. [TJ 
Using properties of functions of several variables 1341 167j . these diagrams can be expressed in 
terms of hypergeometric functions of one variablqlj, whose e expansions up to weight 4 are 
presented in Ref. [56l ISU |66] El and available via the web [70] . We recall that up to weight 
4, the e expansions of all master integrals collected here are expressible in terms of Nielsen 
polylogarithms only. The hypergeometric representations have been derived also in [68j for Ci 
and C2, in [281 [67] for Cq and in [26j for Cn. Up to the finite part, some of these diagrams have 
been studied in [S^]. For certain diagrams (C4, Ce, Cg, Cio, Cn), the e expansion of the first 
several coefficients was given in Ref. |42j in terms of multiple polylogarithms of two variables. 
We use the notations ji23=ji+j2+j3, jmn = jm+jn below. 

We will conclude with a review of the results for special cases: 

• The massless triangle diagram with one massless external on-shell momentum is express- 
ible in terms of two Gauss hypergeometric functions. This result follows directly from a 
relation in Ref. [26]. The Cayley determinant vanishes in this case. 

^"^We are indebted to A. Davydychev for assistance in that analysis. 
^""This is enough for the calculation of two-loop corrections. 
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The analytical result for diagram C\ with arbitrary powers of the propagators is express- 
ible in terms of a Gauss hypergeometric function with one integer upper parameter: 



Ci 



2\n 



/2-,,... r(ii23-t)r(f-ii3) 



n 



J123-2) Jl 

n 
2 



r(f)r(i2) 

The differential reduction will result in one Gauss hypergeometric function. The Cayley 
determinant vanishes for Ci. 

The diagram C2 with arbitrary powers of propagators is expressible in terms of two 
hypergeometric functions 3F2. In this case, both the Gram and Cayley determinants are 
nonzero, and the master integral is 



Co 



ITT' 



n/2 



-(m^)2 



-3jr(3-f)r(f-2) 



r(f) 



2F1 



+ I ' ^'^-}^^^-^> 2F1 



r(n-2) 



1,1 

n 
2 



1 - - 1 

n- 2 



2 



m 



For diagram C3, the result for arbitrary powers of propagators is expressible in terms of 
the function 3^2. Both the Gram and Cayley determinants are nonzero, and the master 
integral is a combination of two Gauss hypergeometric functions: 



C3 

n/2 



= -(m^)2 



n 3r(f-2) rr(n-4) 



r(n-3) 1 r(§-i) 



+ 



Q2X 2-2r(f-i)r(2-f) 



r(3-f) 



2F1 



1,1 

5 — n 
I - - 1 

Q n 
2 



The diagram C4 with arbitrary powers of propagators is expressible in terms of a Gauss 
hypergeometric function with one integer parameter: 

C4 _ r(ii23-f)r(§-ii3)r(n-ii2-2j3) 



lb 

(_^2yi23-2 r (n-ji23) r (§-j3) r(j2) 



2-^1 





2 5 Jl 


91] 


("1 


-J3 


m? J 



For arbitrary powers of propagators, the diagram C5 is expressible in terms of the Appell 
function Fi: 



r(ji23-§)r(f-ji2) 



P ( ^- n ■ A . n Ql Q2 

-''1 Jl23-2 5jl5j2,2 ~2'~2 



il-n^n/2 V - / - r(i3)r(f) 

When the squared external momenta are equal, Qi = Q2 = Q^, it reduces to the Gauss 
hypergeometric function: 

C5 , .__2^g-,•,z■,^(il23- 



^1— nyj-n/2 



(_^2^ 2-^123 



n\ -p In 
2) ^ \2 



-J12) 



r03)r f 



2-f 1 



J123 - 2' J12 

n 
2 



n=Qi=Q^ - V2. 

For Qi = Q2, the Gram determinant is zero, and when Ql = Q2 = m'^, the Cayley 
determinant is also zero. 

For Cg, both the Gram and Cayley determinants arc nonzero, and 



n/2 



ITT 



-(m^)2 



.3[ r(3-f)r(n-5) 
1 r(n-3) 



+ 



Q2^2-2 r2 (n_i)r(2-t; 



r(n-2) 



2-^1 
3 - 



1,1 

7-n 
2 



4m^ / 



n 



1,' 



Aw? ) 
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• The diagram C7 with arbitrary powers of propagators is expressible in terms of the 
function 3^2- For this diagram, both the Gram and Cayley determinants are nonzero, 
and the master integral is 



C7 



n/2 



n 3r(t-i)r(3-f) 



r(i) 



■iF2 



1,1,3 



2 ' ■ 



• The diagram Cg with arbitrary powers of propagators is expressible in terms of the 
function 4F3. For this diagram, both the Gram and Cayley determinants are nonzero. 
The master integral is 



n/2 



[m 



2,^_3r(§-i)r(3-§) 



r(§) 



3^2 



1 o n n 

i, O 

n 3 
2 ' 2 



1 



Am? J 



For Cg, both the Gram and Cayley determinants are nonzero. 



n/2 



(rn^) 2 



i^-3r(3-§)r(§-i) 1 



r(f) 



iF2 



3-§,l,l 

n cy 
2 ' ^ 



Q\ Q2 
3F2 



3-1,1,1 

2'^ 




When = the Gram determinant is equal to zero. 

For diagram Cio, the Cayley determinant vanishes, so that the diagram can be reduced to 
a linear combination of one-loop propagator diagrams (see Ref . ^37j ) . The hypergeometric 
function representation is 

Cio _ r(3-f) 



2Q2(n-4) 

xJ (Q2 + ^2_^2)(^2)§-3 



2^1 



1,3-n 

2 



((52 + ml- 
((52 - + mi)2 



1,3- 



• For this diagram, both the Gram and Cayley determinants are nonzero. The master 
integral is 

/ — - I I i^^2 

3F2 



C 
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1 



ITT 



■n/2 



(m2)§-3r (3- 



3-1,1,1 

2' ^ 



4m 2 



The all-order e expansions of Cn is expressible in terms of multiple polylogarithm of a 
square root of unity. 

The master integral for diagram C12 was evaluated in Ref. [67j in terms of a linear 
combination of two 3-F2 functions of the same type, as for the diagram Cg: 



C 



12 



n 

in 2 



n 



1 



2 J 2{Ql-Ql 



iF2 



3-|,l,l 
2 ' ^ 



Ql] 

Am? ) 



Ql - 3F2 



3-t,l,l 
2 ' 



— ] Q^ 
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